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Abstract 

In this paper we study the general reconstruction of a compactly supported function 
from its Fourier coefficients using compactly supported shearlet systems. We assume that 
only finitely many Fourier samples of the function are accessible and based on this finite 
collection of measurements an approximation is sought in a finite dimensional shearlet re¬ 
construction space. We analyse this sampling and reconstruction process by a recently 
introduced method called generalized sampling. In particular by studying the stable sam¬ 
pling rate of generalized sampling we then show stable recovery of the signal is possible 
using an almost linear rate. Furthermore, we compare the result to the previously obtained 
rates for wavelets. 


1 Introduction 

A general task in sampling theory is the reconstruction of an object from finitely many mea¬ 
surements. Such a problem also appears in more applied fields such as signal processing and 
medical imaging, the latter being a main motivation of this paper with view to applications. A 
Hilbert space where such problems are often modeled in is the Hilbert space L 2 (R 2 ) - the space 
of square integrable functions -, indeed, the reconstruction task in magnetic resonance imaging 
(MRI) is mostly modeled in L 2 (R 2 ). In this paper we assume the object of interest can be 
modeled by a compactly supported function / G L 2 (R 2 ) and the acquired samples are simply 
the Fourier coefficients of this function. Of course, a reconstruction of such a function / can 
then be obtained by forming a truncated Fourier series using the sampled Fourier coefficients. 
However, physical constraints make it impossible to sample all Fourier coefficients which causes 
an error in the reconstruction and possibly leads to non satisfactory reconstructions. For ex¬ 
ample the Gibbs phenomenon is a typical artifact that is strongly visible in the approximation 
if not enough Fourier coefficients are used in a Fourier series. Sampling Fourier coefficients 
and forming a reconstruction by a Fourier series means that sampling and reconstructing the 
signal is both done by using the same system. While the way we acquire the data is indeed 
usually fixed by the acquisition device, for instance in MRI, the reconstruction system is not. 
Very often additional information about the object is known, e.g. the shape or the regularity. 
Therefore it is plausible to use this extra information in the reconstruction scheme by choosing 
a suitable reconstruction system that favors a representation of this object. 
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1.1 Sampling and reconstruction model 

As already mentioned at the beginning of the introduction, many reconstruction problems can 
be modeled in a Hilbert space T~L with inner product (■,•). The samples are then assumed to 
be acquired with respect to a fixed sampling system {si, S 2 , • • •} C H and the reconstruction 
is computed using another system {n,r 2 ,...} C 1~L which can be arbitrarily chosen. The 
concrete reconstruction problem can then be formulated as follows: Given finitely many linear 
measurements G N of a function / G T~L, find a reconstruction /jv G 

spanjri,..., rjy}, N G N that is numerically stable and whose error converges to zero as N 
tends to infinity. 

Reconstruction problems of this type have already been studied, e.g. by Unser and Aldroubi 
in (34| and later also by Eldar in ns mg for the case N = M where the used method is 
called consistent reconstruction. We investigate this reconstruction problem by using generalized 
sampling (GS) which was introduced by Adcock et al. in a series of papers O0li6j. In fact, GS 
can be seen as an extension of consistent reconstruction and a comparison of both methods is 
presented in [6j. The key difference of GS to consistent reconstruction is to allow M and N to 
vary independently of each other. This additional flexibility can be used to overcome some of 
the issues that are addressed in [6]. GS further provides a quantity called the stable sampling 
rate, which represents the number of samples M G N that are needed in order to obtain an 
approximation that is stable and convergent using the first N € reconstruction elements of 
the a priori chosen and ordered reconstruction system. We will introduce generalized sampling 
and in particular the stable sampling rate in sufficient detail in Section [2j 

1.2 Sparsifying systems 

The selection of a different system for reconstruction than the given sampling system can only be 
done faithfully if the structure of the underlying object is known. Indeed, it is well known that 
most natural images, in particular medical images, have a sparse representation with respect 
to a wavelet basis |12j . meaning either many wavelet coefficients are zero or have a fast decay 
with respect to increasing scales. This sparse representation by wavelets is one of the reasons 
why they became very important in imaging problems, see e.g. [33j. 

However, over the years many other sparsifying systems similar to wavelets have arisen such 
as curvelets [9j and shearlets |24| . More importantly, these systems are known to outperform 
wavelets in higher dimensions in the sense that both, curvelets and shearlets provide a so-called 
almost optimal sparse approximation rate of cartoon-like images which are a standard model 
for natural and medical images, see H113I25], while wavelets provable do not reach this rate. 
For this reason, it is of utmost importance to study not only wavelets but in addition shearlets 
as a reconstruction system to approximate a compactly supported function based on its Fourier 
coefficients in the context of generalized sampling. Indeed, wavelets has been studied in mu 
see also Subsection 11.31 

Note that there exists a general notion of parabolic molecules introduced by Grohs et al. 
|19| for which shearlets and curvelets are special instances. A natural question is then, whether 
our results can be extended to such systems and we will comment on this in Section [3 

1.3 Contribution and related work 

To enhance the quality of the Fourier inversion of the sampled data, great effort has been done 
in applications to improve the reconstruction process by changing the reconstruction basis or 
adapting the sampling process. For example in [21] the authors introduced wavelet-encoding 
which allows the user to directly sample wavelet coefficients instead of Fourier coefficients by 
adapting the acquisition device. However, as for now this method also has several shortcomings 
such as low signal-to-noise ratio [30l| . On the other hand, the proposed method GS can be 
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understood as a post processing method that does not require to change the acquisition device 
since the sparsifying transform will be incorporated in the reconstruction process and not in 
the sampling process. 

The use of GS to study the reconstruction problem from Fourier measurements has already 
been investigated for wavelets in ID [7] and 2D [5]. The authors of the respective works 
showed, that for dyadic multiresolution analysis (MRA) wavelets the stable sampling rate is 
linear, meaning up to a constant one has to sample as many Fourier coefficients as one wants 
to reconstruct wavelet coefficients. Moreover, it was shown in [7] that a defiance of this rate 
leads to unstable approximations. 

The main result of our paper shows that the stable sampling rate is almost linear for sam¬ 
pling Fourier coefficients and reconstructing shearlet coefficients. The main challenge compared 
to the previous results of compactly supported MRA wavelets is that shearlets neither obey 
a classical scaling equation which enables an expansion of translations of a single function at 
highest scale nor form an orthonormal bases. However, both properties are heavily used in the 
proofs for the wavelet case. Therefore the proof requires a different strategy and is build upon 
good localization properties in frequency of the reconstruction system. Moreover, to the best 
of the authors knowledge there are so far no results for a stable sampling rate for redundant re¬ 
construction systems that merely form a frame. In fact, the redundancy or instability measured 
by the lower frame bound of the reconstruction system now becomes apparent and effects the 
stable sampling rate directly. Such a penalty does not occur when using orthonormal bases or 
Riesz bases. Therefore the rate that we obtain for shearlets is worse than the rate for wavelets. 

Since the stable sampling rate is worse for shearlets than for wavelets we justify the theoret¬ 
ical gain of using shearlets in combination of GS by showing that the error of the GS shearlet 
reconstruction converges asymptotically faster to zero than the error of the GS wavelet recon¬ 
struction under some assumptions if both reconstructions are computed from the same number 
of measurements. 

As for practical applications such as MRI, wavelets are often used as a reconstructions system 
and it has been shown that these often lead to superior results in terms of lower sampling 
rates and reduced artifacts |28| . We will provide some numerical tests with MR data and 
deterministic subsampling patterns in the numerics section using shearlets and compare these 
to reconstruction methods based on wavelets and Fourier inversion to show the strength of using 
shearlets as a reconstruction system for the recovery problem from Fourier measurements. 

1.4 Notation 

In GS the sampling system and the reconstruction system are allowed to be almost completely 
arbitrary, in fact, assuming the frame property is enough. Recall that a countable sequence 
( fi)iel i s called a frame (fl'lj) for a Hilbert space Ti with inner product (-,-) if there exist 
0 < A < B < oo such that 

A\\f\\ 2 <J2\(fJi)\ 2 < B Wf\\ 2 , V/eH. 

i£l 

For most parts we will assume H to be the space of square integrable functions L 2 (M 2 ) equipped 
with the standard inner product (/, g) := f f~gd[i and induced norm ll/ll : = >/(/, /> for /,se 

l 2 (k 2 ). 

The Fourier transform of a function / G L 1 (]R 2 ) is denoted by / and we use the definition 

f(0= I f(x)e~ 2 ^dx, £ G K 2 , 

J R 2 

and with a standard approximation argument we extend this definition to functions in L 2 (R 2 ). 
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We will also often write to indicate an inequality to hold up to some constant that is 
irrelevant or independent of all important parameters. Similarly "x" should denote equality 
up to some constant. 

1.5 Outline 

We continue with a short review of generalized sampling in Section [2j Then we introduce 
shearlets and recall the concept of sparse approximation of cartoon-like functions in sufficient 
detail in Section [3J Section []] contains the precise definitions of the sampling and reconstruction 
spaces that are used for the main theorem that is formulated in Subsection 14.31 All proofs are 
then presented in Section^ Finally, in SectionOwe aim to numerically confirm the methodology 
and effectiveness of using shearlets as a reconstruction system. 

2 Generalized sampling 

Generalized sampling was introduced by Adcock et al. in mm® as a reconstruction method 
for almost arbitrary sampling and reconstruction systems. In fact, similar to consistent recon¬ 
struction (Hi cna \m these systems are assumed to form a frame. As already mentioned in 
the introduction, the key difference of generalized sampling and consistent reconstruction is 
to allow the number of measurements to vary independently of the number of reconstruction 
elements. This flexibility enables generalized sampling to overcome some of the barriers of 
consistent reconstruction [3]. 

2.1 GS reconstruction method 

Given a sampling system {si, S2, ■ ■ •} C 7~t we define the sampling space S C If as the closure 
of its span, i.e. 


S = spanjsfc : k € N}. 

The finite dimensional version of this sampling space is denoted by 

Sm = spanjsi,..., %}, Me N. 

The sampling vectors {s^ : k € N} are used to model the measurements, in fact, we assume 
that the samples are given as linear measurements of the form 

m f {k) = ( f,s k ), k € N. (2.1) 

Analogously for a reconstruction system {r±,r2, ■ ■ ■} C T~i we define the reconstruction space 
1Z C % as 


TZ = span{rfc : k G N} 

and, likewise, its finite dimensional version is defined as 

TZn = spanjri,..., r^v}, JVeN. 

We will assume the sampling system to form an orthonormal basis and the reconstruction system 
to be a frame, since this will exactly be our setup, see Section [2 A fully general presentation of 
generalized sampling can be found in Si]. Further it is natural to assume a certain subspace 
condition, indeed, it is required that 

1Z fl 5^ = {0} and 1Z + S is closed. (2.2) 
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This guarantees a well posedness of the finite dimensional reconstruction problem, cf. Theorem 
12.11 Indeed, the reconstruction problem is now formulated as follows: given a finite number of 
measurements (/, si),...,(/, sm) of some unknown / G P we wish to determine a reconstruction 
/n € TZn such that ||/ — /jv|| is small and f]\r —> f as N — > oo fast. The following theorem 
leads to the generalized sampling reconstruction and proves its existence. 

Theorem 2.1 ([6]). Let Sm and TZn be as above and Ps M be the following finite rank operator 

Ps M '■ 7-1 Sm 

M 

f ^ E (/j &k} Sfc- 

k =i 

If (12.2(1 holds, then there exists an M £ N such that the system of equations 

{Ps M fN,M,rj) = {Ps M f,rj ), j = 1, • • • ,N (2.3) 

has a unique solution f^,M € TZn- Moreover, the smallest M £ N such that the system is 
uniquely solvable is the least number M € N so that 

c N ,M ■■= inf \\Ps M f\\ > 0- ( 2 - 4 ) 

/ t K'N 5 
11 / 11=1 


Furthermore, 

11/ - Pn N f\\ < 11/ - /at,m|| < —11/ - iWH, (2-5) 

cn,m 

where Pn N : Ti —>• TZn denotes the orthogonal projection onto TZn- 

Definition 2.2 ([S]). The solution }n,m ' m Theorem \2.1\ is called generalized sampling recon¬ 
struction. 

Definition 2.3 ([351133]). The quantity cn,m in (El is called the infimum cosine angle between 
the subspaces TZn an( T Sm- 

2.2 Stable sampling rate 

In light of Theorem 12.11 it is crucial to have control over cn,m- This motivates the definition of 
the stable sampling rate. 

Definition 2.4. For any fixed N € N and 6 > 1 the stable sampling rate 0(iV, 6) is defined as 


The stable sampling rate determines the number of measurements M that are needed in 
order to find an approximation /jv,m € TZ n or equivalently, find N coefficients such that the 
angle cn.m is controlled by the threshold 9. We wish to emphasize that the control of cn,m 
allows the control of the error that is made by the approximation via GS by (12.5|) . 


0(iV, 9) = min < M € N : cn,m > 
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3 Shearlets 


In this section we introduce the function systems that are used for the reconstruction space, 
these are compactly supported shearlets. 

Shearlet systems were first introduced by K. Guo, G. Kutyniok, D. Labate, W.-Q Lim and 
G. Weiss in EDI m and the following notations have became standard in this topic. The 
parabolic scaling matrices with respect to scale j € N U {0} are denoted by 

— f 2j 0 \ 7 / 2- ? / 2 0\ 

2j Vo y/ 2 ) 1 2j \ o y) ’ 

and the shearing matrices with parameter k £ Z are 



These operations, together with the standard integer shift of functions in L 2 (M 2 ) are then used 
to define the cone adapted discrete shearlet system. 

Definition 3.1 ([22]). Let q £ T 2 (M 2 ) be the generating functions and c = (ci,C 2 ) £ 
(M + ) 2 . Then the (cone adapted discrete) shearlet system is defined as 

STL= $(<)>, ci) Uf(<),c) U$(^,c), 

where 

$(</>, ci) = {</>(■ — c\m) : m € Z 2 }, 

and 

T(i/gc) = {Vh,fc,m = 2 3j/4 i/> ((Sfc^aj) • -cm) : j > 0, |fc| < 2 i/2 ,m € Z 2 } , 

W,c) = {Vh'.fc,m = 2 3j/4 '0 ((SIA 2j ) ■ —cm'j : j > 0, \k\ < 2 j/2 ,m£ Z 2 } , 

where the multiplication of c and c with the translation parameter m should be understood entry 
wise. 

The shearlet system has a multiscale structure, however, it does not form a classical mul¬ 
tiresolution analysis (MRA), see |12l!29j for a definition of an MRA. Moreover, it is still an open 
question whether there exists compactly supported orthonormal shearlet bases. Nevertheless, 
shearlets can form a frame. 

Theorem 3.2 ([22]). Let £ 2 (M 2 ) such that 

|<K6,6)1 < Ci min{l, |£i| -r }min{l, |6l~ r } 

and 

1^(6,6)1 < C' 2 min{l,|:6| a }min{l, |6r r }min{l, |6r r }, 

for some constants C±,C 2 > 0 and a > r > 3. Further let ip(xi,X 2 ) = if(x 2 ,%i) an d assume 
there exists a positive constant A > 0 such that 

ife)i 2 + E E lAsAtr'oi 2 + E E i5(s*A-)- 1 oi 2 > a 

J>0 |fe|<[2J/2] j>0 |fe|<[2J/ 2 l 

holds almost everywhere. Then there exists c = (ci,C 2 ) £ (M + ) 2 such that the cone-adapted 
shearlet system forms a frame for L 2 (M 2 ). 

An explicit construction of compactly supported shearlets that form a frame can be found 
in [22] • Note that the frequency decay assumptions on the generators given in Theorem 13.21 are 
essential and secure that the overlapping in frequency domain is controlled. In particular, the 
frame constants depend on this behavior, [23 GUI!!]. For more information about shearlets we 
refer to [23 ] and the references therein. 
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3.1 Sparse approximation of cartoon-like functions 

Shearlets are vastly used in applied harmonic analysis due to the optimal sparse approximation 
rate for cartoon-like functions. We now recall the definition of cartoon-like functions and the 
results about optimal sparse approximation rate. 

The set of cartoon-like functions was first introduced in |l3| and since then has been fre¬ 
quently used as a mathematical model for natural images. 

Definition 3.3. Let v > 0 and f : R 2 —> C a function of the form 


f = 9 + hxB, 


where B C [0, l] 2 is a set whose boundary is a closed C 2 curve with curvature bounded by v 
and g,h E C 2 (R) are compactly supported in [0, l] 2 with H^Hc 2 ) H^-Hc 2 < 1- Then we call f a 
cartoon-like function. The set of cartoon-like functions is called class of cartoon-like functions 
and is denoted by £ 2 (n). 


The following theorem shows the sparse approximation rate for shearlets. 
Theorem 3.4 1 [25] ). Let S'H(4>,if,if,c) be a compactly supported shearlet frame with 

i) |-0(C)| < C'min{l,|^i|“}min{l,|^i|- r }min{l,|6r r }, 


a) 


9 m 




^ IM£i)l (i + 


where a > 5, r > 4,h E T 1 (R), C > 0 is some constant. Further assume if fulfills i) and ii) with 
switched roles of f \ and £ 2 - Then for any fixed v > 0 and f E £ 2 {v) we have 


||/-/;v|||<C'(logfV) 3 iV- 2 , as N —> 00 , 


where /at is the best N-term approximation obtained using the N largest shearlets coefficients 
of f in magnitude allowing polynomial depth search and C' is some constant independent of N. 

Remark 3.5. The rate obtained in Theorem \3.f \ is also achieved in /?/ for band-limited shearlet 
systems and for curvelets in W- 

In [T3 j it was shown, that under some assumptions on the representation system and the 
selection procedure of the largest coefficients in magnitude, the optimal approximation rate of 
a cartoon-like function / obeys 


\\f-fN\\ 2 <CN~ 2 . 

This rate is, up to a logarithmic factor, achieved by shearlets, cf. Theorem 13.41 Moreover, 
compared to the N~ 2 factor, the logarithmic term can be neglected and, thus, the rate in 
Theorem 13.41 is referred to as almost optimal sparse approximation rate. 

Remark 3.6. It is commonly known, that standard wavelet orthonormal bases only obtain a 
rate of order N^ 1 , see M- Moreover both rates for shearlets and wavelets are sharp, that is, 
there exists cartoon-like functions such that the rate is achieved. 

3.2 Assumptions on the generator 

For the rest of this paper we assume (j> and if to be compactly supported functions in L 2 (R 2 ) and 
define if(x 1 , 2 : 2 ) := ' l f( x 2 ,xi)- Moreover, we assume the shearlets to have sufficient vanishing 
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moments and frequency decay, more precisely, we assume there exist some constants C \, C 2 > 0 
such that 

&)l < C 1 ' (l + |^|)r (l + |^ 2 |)r’ 

and 


l<A(6 


1^(6,6)I < C 2 ■min{l,|fi| a } ■ 


1 1 

(1 + I 6 l) r (1 + l£ 2 |) r ' 


(3.2) 


where the regularity parameters a, r > 0 are large enough so that the shearlet system forms a 
frame for L 2 (R 2 ). These assumptions are equivalent to the frequency decay assumptions stated 
in Theorem 13.21 


4 Stable shearlet reconstructions from Fourier measurements 

We will now define the reconstruction space and the sampling space for which we determine a 
stable sampling rate. 

4.1 Shearlet reconstruction space 

Without loss of generality we can assume that the generating scaling function (f> and the shearlets 
ip, ip are compactly supported in [0, a ] 2 where a is some positive integer. Then we consider all 
scaling functions whose support intersect [ 0 , a ] 2 and denote this index set by 17, i.e. 

17 = {m € Z 2 : supp cp m n [0, a ] 2 / 0} = {(mi, m 2 ) G Z 2 : —a < mi, m 2 < a}. 

Similarly, we consider all shearlets whose support intersect [0, a] 2 . For this, let J — 1 G N U {0} 

be a fixed scale. Then A j denotes the following paramater set 

Aj = | (j. k, m ) G Z x Z x Z 2 : 0 < j < J — 1, |fc| < 2 J / 2 , m G ^7j,fc| , 

where 17 j k = {m G Z 2 : supp ipj,k,m FI [0, a .] 2 7 ^ 0} is, due to the compact support of each 

(• tpj t k,m)m j °f finite cardinality and for i/j and the second cone we analogously write 

Aj = |(j) k, m) G Z x Z x Z 2 : 0 < j < J — 1, \k\ < 2 ^ 2 ,fh G 

with !7j£ = {m G Z 2 : supp F [0,a ] 2 7 ^ 0} being of finite cardinality, respectively. In 

fact, the scale determines the number of translation and will be of order 2 3 / 2j for a fixed scale 
j. The reconstruction space at scale J — 1 is then defined as 

ll Nj := span| {cj) m : m G 17} U {ipj,k,m ■ ( j,k,m ) G A,/} U : ( j,k,m) G Ajj 

(4.1) 

Up to a fixed scale J — 1 we, asymptotically, have Nj = 2 2J many generating functions in 
TZnj as J —> 00 . Note that by construction at each scale we only have finitely many elements, 
therefore, an ordering can be performed quite naturally namely we order the system along scales 
and within the scales, the shearing is ordered from — 2 J / 2 to 2 J//2 and the translations in Z 2 are 
ordered in a lexicographical manner for each fixed shearing and scaling parameter. 

By TZ we denote the reconstruction space that contains all shearlets across all scales, i.e. 

| {4> m : m G 17} U {ip jtk ,m ■ ( j , k, m) G Aj, J G N U {0}} 

U : e Aj, J G NU {0}| j. 
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1Z = span 





4.2 Fourier sampling space 

To define the sampling space we first choose T\ , T2 > 0 sufficiently large such that 

TZ C L 2 ([—Ti,T 2 ] 2 ). 

Let £ < t 1 +t 2 < ^ con t r °l the sampling density. Then we define the sampling vectors on the 
uniform grid Z 2 by 

=^ 2 ^->-X[-t 1 ,t 2]2j ^€Z 2 . (4.2) 

Based on these sampling vectors we define the sampling space S by 

= span |4 £) : l € Z 2 j . 


For M = (Mi, M2) S N x N let 

4? = span{4 £) : l = (h,h) £ T?,-M { <ti<M i ,i = 1,2} 

be the finite dimensional sampling space. Note that M = (Mi, M2) determines the size of the 
measured grid and the total number of possible measurements are in this case asymptotically 
of order Mi • M2. 

The task for us is now to describe the relationship between M and N in terms of the stable 
sampling rate which means such that stable and convergent reconstructions exist. In the event 
that the stable sampling rate is linear, i.e. ©(IV, 0) = 0(N) we would have Mi • M2 = 0(2 2J ) 
as for the wavelet case, mm- 


4.3 Main results 

Our main result shows that the angle between the shearlet reconstruction space and the Fourier 
sampling space can be controlled with an almost linear stable sampling rate. 

Theorem 4.1. Let STL^jif,^) be a compactly supported shearlet frame with generators 
and if and let N < Nj = 0(2 2J ). Then for all 9 > 1 there exists Sq > 0 such that 

c N , M = inf \\P s (e)f\\>L 

f$zR,N ^ M U 

!l/ll=i 


where M = (Mi, M2) € N x N with Mi = \SqA^^ 2 ' ^ • 2' 7 ^ 1+ ^/e], 8 > an d r > 0 is the 

regularity parameter from \3.1\ ) and H3.2\) . Therefore Q(N,6) = 0(N 1+s A N L ^ 2r ^). Further, 
the constant Sq does not depend on N but on 9, a, and r. 


Remark 4.2. The rate appearing in Theorem \4-l\ is slightly worse than the rate for wavelets 
where the stable sampling rate is indeed linear, i.e. Q(N,9) = O(N). More precisely, one 
has to oversample by a factor IV 5 and the rate depends on the lower frame bound An of the 
finite shearlet system that is used to span TZn- 


Remark 4.3. If we interpret the reciprocal of the lower frame bound as a measure for redun¬ 
dancy or stability, then our result states that more samples are necessary if the system is more 
redundant or less stable. Note that if the frame is tight, then the additional factor Af^ 1 can be 
neglected. 
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If one uses the same number of shearlets and wavelets to reconstruct a function from its 
Fourier measurements, then by Theorem 14.11 one has to acquire more Fourier measurements 
to guarantee the existence of the GS solution if shearlets are used compared to wavelets. This 
immediately raises the question why the use of shearlets in the context of GS is justified, if 
one has to acquire more samples to guarantee stable and convergent solutions. We address 
this issue from a different perspective. Suppose we were to sample an image using the rate 
computed in Theorem 14.11 for shearlets, but we want to reconstruct the image from the same 
number of measurements using shearlets and wavelets. Then, due to the stable sampling rate 
for shearlets and wavelets, we are allowed to use more wavelets than shearlets. More precisely, 
let a : N —> N be the following oversampling function 


a(N) 


N 


1+5 


1 


2 

A 2r— i 

A N 


N £ N. 


(4.3) 


Then the following result holds. 

Corollary 4.4. Let be a compactly supported shearlet frame with generators (f>,if, 

and if with sufficiently large regularity r and 6 > 2/{2r — 1). Further, let f be a cartoon-like 
function f. For JV G N denote by fff the best N-term approximation of f using shearlets 
and f™ {N) the best a (TV)- approximation using wavelets, cf. Subsection 1 3.11 Further, let Nj = 
2 2J , J £ N be the smallest number such that N,a(N) < Nj. If 1 Z s Nj denotes the shearlet 
reconstruction space space so that fff £ F s Nj and denotes the wavelet reconstruction 

space so that f™^ £ 7 then 


11/ — G s NjtM (f)\\ < TV -1 (log TV) 3 / 2 , 


and 


\\f-G w a{Nj ), M '{f)\\<< N )- l/ \ 

where M and M' are of order a(Nj ) and G s Nj M (f) and M (f) are the GS solution with 
respect to the spaces F s Nj and LVf , Nj ). 

Remark 4.5. Note that IV -1 (log N) 3 / 2 < cr(A r )" 1 / 2 if 

N - (1 “ 5) / 2 (log TV) 3 / 2 < A]i {2r - 1] , (4.4) 

and therefore, if ||/ - /^|| x \\f - G s Nj M {f)\\ and \\f - /^|| x \\f - G^ Nj)M {f)\\, then 

11 / — Gr s NjM (f)\\ < ||/ — Gff^ Nj)M ,{f)\\, 

where "x " means equality up to some constant. 

In general, there is no control about the lower frame bound of the finite shearlet system 
known, thus, the reader might wonder whether (14.41) can hold. In Figure [2] in Section [U] we 
numerically confirm the validity of such an behavior of the lower frame bound. 


5 Proofs 

We first sketch the idea of the proof of Theorem 14.11 
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5.1 Intuitive argument for the proof of Theorem 14.11 

In order to prove the main theorem we have to bound 


inf 

f£K N 

11 / 11=1 



from below for respective M and N. For any / £ TZn we have 

ii^/ii 2 = E i</.4 e| )i 2 , 

where Im = = (^ 1 ,^ 2 ) £ Z 2 : —Mi < £i < Mi,i = 1,2}. Hence, we could equivalently bound 

l|ri.,/ll 2 = E K/.ofil 2 , P-1) 

from above. Now, the key idea is to use the effective frequency support of the shearlets, in 
particular, the energy of ( 4’j,k,m) A is essentially localized in frequency bands of width 2 2j up to 
some constant, see Subsection 15.21 for precise statements. Moreover, since f € ~R-n is a linear 
combination of shearlets up to a fixed scale J — 1 the function / might be essentially supported 
in [— C2 J , C2 J ] 2 for some constant C as well. Hence, by making the grid Im large and dense 
enough the term in (15.ID should be small. However, there are two main concerns: First, the 
effective support could grow faster than linearly with respect to scaling factor 2 3 . Second, 
taking linear combinations of shearlets might destroy the control of the frequency behavior and 
hence the essential supports in Fourier domain. Unfortunately, both obstacles are inevitable 
causing the additional factor N s and the dependence on the lower frame bound An in the stable 
sampling rate, respectively. More precisely, if / £ IZn then there exists ci,..., cjv £ C such 
that 

./' = X/ Cx ^ x 

A <N 

but in general we have no control over the coefficients (c\)\< n- Therefore in order to resolve 
the latter concern, we expand / as follows: 

/ = ^2 if, S^Va^a- 

A <N 


Then 


K/» 5 jvV a) I 2 

A <N 



This is how the lower frame bound of the finite shearlet system comes into the proof and thus 
effects the stable sampling rate. Note that this problem does not appear if the reconstruction 
system is an orthonormal basis or a Riesz basis, thus, an additional penalty of the stable 
sampling rate is not necessary for wavelets. 


5.2 Effective frequency support 

The decay assumptions (13.11) and (13.21) give rise to the essential support of each shearlet atom 
ip\ in frequency, see Figure Q] The next result secures that the effective frequency support is 
controllable. 
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Proposition 5.1. Let J £ N and (ipx)x<Nj be all shearlets up to scale J — 1 and uj > 0. Then 
there exists a constant S := S(uj,r,e) such that for Im := {(^1,^2) € Z 2 : —Mi < £i < Mj, i = 
1,2} with Mi = S2 J ( l+5 \i = 1,2 and 5 > we have 

E E l(*0 A (rf)l 2 £ 

te(/ M ) c A<Afj 



Figure 1: Effective frequency support of shearlets and tiling of the frequency plane, light regions 
correspond to large values, therefore regions with large energy and dark regions to small values, 
i.e. regions with small energy. Images are computed using the ShearLab package downloaded 
from http: //www. shearlab. org/ 


To prove Proposition 15.11 we will use the following lemma which can be found in |17| (Ap¬ 
pendix K). 


Lemma 5.2 (EH). For y € R, r > 1, and a,b > 0 we have 

1 1 


r°° 1 

Jo (1 + a 


■dx < 


1 


a|x|) r (1 + b\x — y\) r ~ max(a, b ) (1 + min(a, b)\y\) r ' 

Proof of Proposition 15.11 Let lo > 0. Then we have to show the existence of an M = (Mi, M2) € 
N x N, so that 


E E + E E + E E E i(« A (rf)i 2 < 

Jb|>M 2 |5|>Mi \h\<M 2 \l 1 \>M 1 \h\>M 2 \h\<M 1 ) \<Nj 


(5.2) 


where Mi, M2 scales like l S , 2 J ( 1+<5 ) with a constant S independent of J sufficiently large. Note 
that 


E i(« A (rf)i L 

\<Nj 


| ee - 2 ^’ m,> ^(d) + 

m'^Q ( j,k,m)£Aj 


235/4 


e -2 *ie((S k A 2j ) T l,m)^ h ( S k A^)~ T z) 


+ 




C j,k,m)eAj 
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2 


2 J_1 27/2 2 J_1 2j/2 ~ _ _ T 

<c*o| + E E |^(Es^r T 0| +E E *•U (vw) 1 

j=0 k=-2i/ 2 j=0 fc=_2i/ 2 V 


since the cardinality of Qyfc and is of order 2 3 / 2j . 
Denote by I, II, III the following terms 


I; - E E + E E + E El 

\|Z 2 |>M2 |Zi|>Mi |Z 2 |<M 2 |Zi|>Mi |Z 2 |>M 2 |Zi|<Mi/ 

/ \ J-l 2^/ 2 

ic= E E + E E + E EEE 

\|Z2|>M 2 |Zi|>Mi |Z 2 |<M 2 |Zi|>Mi \l 2 \>M 2 \l 1 \<M 1 J j=0 k=-2i/ 2 

( \ J -1 2?/ 2 

nc= E E + E E + E E EE 

\\l2\>M 2 \li\>M 1 \1 2 \<M 2 |Zi|>Mi |Z 2 |>M 2 |Zi|<Mi/ J=Q fc=_ 2 7 /2 


(5.3) 


^ ( £ ( ) 1 


-T 


We will estimate each sum by using the decay conditions (13.11) and (13.21) . respectively, in order 
to obtain (15.2p for sufficiently large S independent on J with Mi = 2 — 1 5, i = 1,2. Since 


s4>(el) 


< C\e 2 


(1 + \eh\Y (1 + \d 2 \y 


we immediately have 


E E I ^ ^ E E 


|Z 2 |>M 2 |Zi|<Mi 


< 


|Z 2 |>M 2 |Zi|<Mi 


c? E 


1 (i + M 1 |) 2r (i + M2l) 2 '- 


|Zi|>Mi 


< cf- 


(1 + Mi|) 2r 

1 


(1 + s ,J ( 1 + 5 )) 2r - 1 ’ 

where the constant Ci changed in each step. Analogously, one obtains 

1 


< C 'f 


E E |^ 

\h\<M 2 |Zi|>Mi 

E E |^( d ) 

|Z 2 |>M 2 |Zi|>Mi 


( 1 + S'2 J ( 1 +<5))2r- 1 ’ 


1 


(1 + S2 J ( 1 + s )) 2r 


- 1 ' 


Hence, combinging (15.5[) and (15.61) gives 

I < cf 

Regarding II we first have by (13.2p 


(1 + S’2 J 'H+'5))2r 


-1 


J -1 2J/ 2 J -1 2 H 2 

E E <E E 


j=0 k=-2i/ 2 


c .? 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


J=0 k=-2i/ 2 


(1 + |e2-JZ 1 |) 2 »-(l + | - ek2~H 1 + e 2-.b%|)2r ’ 


(5.8) 
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and, likewise for III, 


J -1 2 i/ 2 

E E 

j =0 k=-2i/ 2 


$ ( £ ( S l A 2J 


-T 


2 J -1 2^ 2 

sE E 

3=0k=-2H 2 


Cl 


(1 + |e2-^ 2 |) 2r (l + |e2-i/ 2 ^ - ek2~H 2 \) 2r 


(5.9) 


We only continue to estimate (15.8p since an estimate for (15.91) is then obtained analogously. 

We consider two cases. The first case concerns shearlets, that are wavelet-like. More pre¬ 
cisely, we consider shearlets with no shearing first, so parabolically scaled wavelets. 

Case I: Let 0 < j < J — l,k = 0. By direct computations similar to (15.51) we have 


E E Ec 

\h\>M 2 |L|<Mi j =0 
J -1 


1 


1 


^E^ 2 E 


(1 + \£2~Hi\) 2r (1 + \e2~^ 2 l 2 \) 2r 

1 X - 1 


E 




J -1 


2 J 


sEgIt E 


1 


>-o e G + \^ i/2 W 


1 


< vV 2 ?!!_ 

- £ 2 q S2 J'(l+5)-j/2\2r-l 

3 =0 V 


< Co 2 


2 3 /2J 


1 


- 2 £ 2 (l + S’2 J ( 1 /2+5))2r-l' 

where C 2 changed over time. In the same manner we obtain 


(5.10) 


J-i 


E E E^ 

N<M 2 |Zi|>Mi j =o 


1 


1 


< Co 2 


, 2 3 / 2J 


(1 + |£2--?7i|) 2t ' (1 + |e2-i/ 2 Z 2 |) 2r _ 2 e 2 (1 + S2 J6 ) 2r ~ l ' 


(5.11) 


and 


J-l 


E E Eg 

|Zi|>Mi |Z 2 |>M 2 3=0 


(l + |£2-i/!|) 2r (l + |e2-i/ 2 ^ 2 |) 2j - 


< C, 


i 2 3 /2J 


1 


2 e 2 ((1 + S2 JS )(1 + S2 J W 2 + s ))) 
Case 2: Let 0 < j < J — l,k ^ 0. By using Lemma 15.21 we have 

J-i 2 H 2 


2i—1 • 


E E E E G 

\h\>M 2 \h\<M 1 j=0 fc = _ 2 j/a 


1 


1 


(1 + |e2-i/i|) 2f - (1 + | - ek2-H 1 + e2~i/ 2 l 2 \) 2r 


1 


J-l 2H 2 

-EE E C 2 4 fc 2 '-i|(l + k2-G||2G/ 2 I 2 /fc|) 2) ' 

j=0 fc= _ 2 j/2 |Z 2 |>M 2 1 IV 1 11 2/12 

k+0 

J-l 2^/ 2 


< 


E E gi 

J=0 k =- 2?/ 2 
fc^O 


2 J |jfe|2J'/ 2 


1 


r|fc| e (l + |^2 J, ( 1 +‘ 5 )-J/ 2 /fc|) 2r - 1 


(5.12) 
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J ~1 2 j/2 9 3/2 j 

E E <£* 


1 


£ 2 n + \S 2 J(l+S)-j/2/ k ^2r-l 

J=0 fc=_ 2 j/ 2 V le ' 

k+ 0 

2 2 2J 1 


£■2 (1 + || 2 J<5 |) 2r - 1 


(5.13) 


Since 


E 


< 


dx < 


2^/2 


Il|<m 2 

we have that 


(1 + | - sk2~Hi + £2-i/2/ 2 |) 2r ~ 7 r (1 + | - ek2-H l + £2~j/ 2 x\) 2r ~ £ 


1 


1 


J-l 23/ 2 

^ ^ ^. ( 62 (1 + |£2 J’Z 1 |) 2r (1 + i - £fc2-J/i + £2-i/ 2 Z 2 |) 2r 

|«2|<M 2 |Zi|>Mi J—0 k=—2 ?/ 2 

fc^O 

J-l 2^/ 2 


< 


EEc 


, 2 3 / 2 - 7 


1 


5=0 k=-Z >/ 2 

2 2 2J 


2 £ 2 _|_ |S2-7(l+tf)-5|)2r-l 


£ 2 (1 + ||2J«|)2r-r 

Finally, the last sum can be bounded as in (15.1311 by first using Lemma [5.21 again 

1 


(5.14) 


1 


J-l 2^/2 

^ ^ ^ ^ (1 + \e2-Hi\) 2r (1 + I -ek2-H 1 + £2-i/ 2 l 2 \) 2r 

|i 2 |>M 2 \h\>Mi J—O / c= „ 2 j/2 
k +0 

J-l 2^/2 „ 3/2l - 1 


< 


E E 


j=o *,_*/* e2 (l + |S»2 J(I+i W|) 2 ’- 1 

fc^O 
o2 J 

< r 2 -_ 

S °2 c2 M , |S( 


1 


£ 2 (1 + || 2 ^|) 2 r-l' 
hence, combining (15.1011 . (15.111) . (|5.12[) . (15.131) (15. 1411 . (15.1511 yields 

1 


o2 J 

n < c|—2 


£ 2 (1 + | f 2 J5 |) 


5oJ(5|\2r-l 


(5.15) 


(5.16) 


Regarding III, we obtain a similar estimate by performing the same computations as for II, 
therefore 

(5.17) 


,2 2J 1 

111 - ° 2 (1 + || 2 J5|)2r-l' 


£ 


,2r—3 


Combining (15.711 . (15.161) . (15.1711 yields 

I + II + III < c'^ 2 j (5 _ 2/(2r _ 1) ^ 2r ._ 1 > 

with a constant that does not depend on J. Therefore, if 5 > 2 r-i the resu lt follows for S 
greater than C'~ 1 ( £2 f 3 ) 1 /( 2r ~i). □ 
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5.3 Proof of Theorem 14.11 

The proof of Theorem 14.11 is a consequence of Proposition 15.11 
Proof of Theorem \f. 1\ Let 6 > 1. Then we want to show 


inf 

11/11=1 



(5.18) 


for appropriate M. For this, let / £ TZn with ||/|| = 1. We prove 


ip-L ri|2 < ^_I 

11 - e 2 5 


(5.19) 


for the claimed M. Since (si)i is an orthonormal system, we have 

ll^(e)/|| 2 = E K/’ S ! £) >| 2 ’ 

where ( Im) c denotes the set complement of Im in Z 2 , in particular 


E k/,4'’>i 2 = E E k/.T’>i 2 + E E i(/.»i‘ ) >i 2 + E E km w >i 2 

\h\>M 2 \l 1 \>M 1 \l 2 \<M 2 \h\>M 1 \h\>M 2 \h\<M! 

Therefore we have 

./ = E VtOV’A = E + E Pj,k,m 1 l } j,k,m + E 


\<Nj 


m'£Q 


(j,k,m)£Aj 


(j,k,m)eAj 


with 


C^m' — (0m'; S/)) Pj,k,m — {f , j,k,mi Si), r Yj,k,m — (&j,k,m’ 1 


Furthermore, 


E 


l a m , | 2 + \Pj,k,m \ 2 + ^2 

m ' en 0,fc,m)eAj (J,fc,m)eAj 


1 '. 


j,k,m 


E K^7,V’a)I 2 <-E 

a N 


\<Nj 


Therefore by Cauchy Schwarz 

E !(/. ®i)l 2 


E ( E + E Pj,k,m' l l J j,k,m+ E 

1£(ImY m ' e ° (j,k,m)eAj (j,k,m)GAj 

E E OL ™ ,( l ) m'{ £ 0+ E 7,fc,mV ; i j',fc,m(e0 + E 

1&(ImY rn'&n ( j,k,m)eAj (J,fc,m)eAj 


< 


E 

1&{ImY 


IN 


^ j 10m' (^0 T ^ j T ^ j £( Pm f 

- m'eO (j,k,m)eAj (j,k,fh)£Aj 


By Proposition 15. II there exists S(6,r,s) so that we can choose 


Mi = 


S(9, r, e)2 J{ - l+ ^ A 2 ^ 1 


£ N, i = 1,2 
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to conclude 


1€{Im) c 

which shows (15.191) and hence 

\\ p s e M f\\ ^ \- 

□ 


5.4 Proof of Corollary 14.41 

For the proof of Corollary 14.41 we recall the following results. By Theorem 13.41 the best AC-term 
approximation for shearlets of cartoon-like functions obeys 

Wf-m <N-\logN) 3 / 2 , 

and for wavelets the best AC-term approximation obeys 

Wf-fN\\<N~ l /\ 

see [23]. Furthermore, these bounds are sharp, i.e. there exist cartoon-like functions such that 
the above inequalities are becoming equalities up to a constant. 

Now, let N € N, fff, P 'a(N J ) he as c n the corollary. By Theorem 14.11 there exist 

for fixed 9 > 1 an M = (Adi, M 2 ) and M' = (M[, M 2 ) in N x N such that 

c Nj ,M = mf ||P S ( £ )/|| > \ 

11/11 = 1 

and 

11 / 11=1 


Moreover, due to Theorem IQ and the linear stable sampling rate for wavelets proven in mm 
we can choose M\ ■ M 2 and M[ ■ M 2 to be of order a(Nj). 

Using Theorem 12.11 we obtain 

11/ - G w a{Nj )Mf )II £ 11/ - P ^ Nj) U )II < 11/ - /“(AO II < 

Analogously, for shearlets we obtain using Theorem 12.11 


2 


\\f-G s NjM U)\\ 2 <\\f 


p n°(m 


< 

rs_/ 


/- (/’^a)iAa 

Ae/jv 


< E k/,</1>i 2 , 

A 0N 


where (V'a)a is the dual shear let system of (i^>\)\ and In denotes the index set of the N largest 
coefficients ((f,ip\)\) with respect to the best AC-term approximation /^-. 

Furthermore, if the regularity is sufficiently large, it was shown in |18| that we can relate 
the analysis coefficients of the primal frame with those of the dual frame by the relation 


where G t is the pseudo inverse of the Gramian operator associated to the shearlet system (i/j\). 
Therefore the claimed rate follows from the decay rate of the shearlet coefficients. Hence, 


\\f-G s NjM (f)\\ 2 < ]T |(/,^ A )| 2 <AC- 2 (logAC) 3 
Hh v 


which yields the result. 
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6 Numerics 


6.1 Inequality (14.411 

In light of Corollary 14.41 we want to compare the bounds IV -1 (log IV) 3 / 2 and a(N)~ 1 / 2 with 
g defined as in (14.311 for N £ N. We therefore aim to verify (14.4[1 by showing the following 
behavior of the lower frame bound An 

Ar -(l- < 5)/2( lo g JV)3/2 < A V( 2r-l)_ (6.1) 

For convenience we will check m in levels, i.e. the lower frame bounds An will always belong 
to full shearlet systems at a scale J. Further, we use the SHEARLAB package which is available 
at 


http://www.shearlab.org 

to generate shearlet systems with respect to different levels and their respective lower frame 
bounds. We wish to mention that both objects are available in the implementation. 

Now, note that a shearlet system consists of approximately many elements for a fixed 
scale J where we ignore constants that are independent of J. We will numerically test (14.4p 
by plotting the lower frame bounds of a full shearlet system for increasing J. This yields a 
monotonically decreasing curve as the scale increase. This will now be compared to the decay 
of the function f(N) = N~^~ s ^ 2 (log IV ) 3 / 2 with N = 2 2J for different values of r and 5 and 
for increasing J. 

The informed reader might notice, that choosing N = 2 2J is significantly moderate, since 
the constant which determines the number of shearlets at scale J is in practice much larger 
than 1. However, for larger N the function / decreases even faster. Moreover, since we are 
only interested in an asymptotic behavior of the lower frame bound with respect to increasing 
scales this is not a restriction. Indeed, already for that case it is strongly indicated that (14.4[) 
hold asymptotically as depicted in Figure [2j 

6.2 Reconstruction from MR data using shearlets 

We further provide some reconstructions from MR data using complex exponentials (Fourier 
inversion), compactly supported Daubechies wavelets and compactly supported shearlets. Al¬ 
though our main result guarantees stable and convergent reconstructions when the sampling 
rate is almost linear it is not efficient (and also not necessary) to acquire that many samples 
in practice. In our numerics we will subsample the Fourier data to achieve practical relevance. 
All computation are done in MATLAB. 

6.3 Data 

The underlying Fourier data or also called k-space was acquired using a multi-channel acquisition 
consisting of four channels. Each of the four k-spaces has a 128 x 128 image resolution. Moreover, 
each channel data results in a single image of same size, e.g. by applying an inverse Fourier 
transform, see Figure [3a] and Figure l3bl 

The single channel images are combined by using the sum-of-squares method |31| . The 
resulting sum-of-squares image from all four channels of the original k-space is used as the 
reference image, cf. Figure l3cl 
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(a) r = 3 and 5 = 0.4,0.5, 0.6 (b) r = 4 and 5 = 0.3,0.4, 0.5 




Scale 


Scale 


(c) r = 5 and 5 = 0.25,0.35,0.45 


(d) r = 6 and S = 0.2,0.3, 0.4 


Figure 2: Plot of the behavior of lower frame bound and the critical function iV l 1 ^^(log JV) 3 / 2 
for different values of r and 5 
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(c) Combined image by sum-of-squares 
Figure 3: Reference images computed with complete data 


6.4 Methods 

We mimic the acceleration of factor 5 by using a 0 — 1 mask to subsample the k-space, cf. Figure 
[4a] and Figure l4bl We have chosen a phy Hot axis spiral which is also used in |35j and a standard 
radial subsampling scheme. 

The final reconstruction is obtained by computing each channel image separately and then 
combine it by using the sum-of-squares method ED- More precisely, if Ik denotes the k -tli image 
reconstructed from the k -tli k-space then the final image I is computed by 

1 = V|Ii| + ... + |/ivi 2 , 


where N is the total number of channels and the the square operation, as well as the modulus 
is to be understood entry wise for the matrices Ik- 

The reconstruction of the single channel images are computed as follows. 

• Shearlets: 

The reconstructions are obtained by solving the following unconstrained minimization 
problem 

min ^||JFu - b\\l + A||5'H(u)||i, 

u Z 

where T is the subsampled Fourier operator, b is the subsampled input data, SH is the 
shearlet transform and A are some parameters larger than zero, and u is the image. 
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The penalty function ||<S'%(it)||i enfores sparsity with respect to the shearlet system. The 
code for the shearlet transform is downloadable at 

http://www.shearlab.org 


• Wavelets: 

Reconstructions are computed using the SparseMRI package from Lustig et al. [28] which 
is downloaded from 

http://www.eecs.berkeley.edu/~mlustig/Software.html 

This method is a standard tool in MRI to compute wavelet sparsity-based reconstructions 
from Fourier measurements. 

• Fourier: 

The reconstructions using complex exponentials are computed by taking the Fourier in¬ 
version of the subsampled data. 

6.5 Results 

The reconstructions show less artifacts and improved image quality using shearlets compared 
to the reconstruction obtained from compactly supported wavelets, cf. Figure [4] Although the 
theory of optimal sparse approximation rates hold in a continuous setup, it is surprising to see 
such strong differences for images of such low resolutions such as 128 x 128. 



Shearlets 

Wavelets 

Fourier inversion 

Spiral mask 

0.638 

0.1006 

0.2570 

Radial mask 

0.0847 

0.1282 

0.2690 


Table 1: Relative error for Figure [4] 
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(a) Phyllotaxis spiral mask: 20.37 % 



(c) Reconstruction: shearlets, spiral 



(e) Reconstruction: wavelets, spiral 



(g) Reconstruction: Fourier inversion, spiral 



(b) Radial mask: 20.74% 



(d) Reconstruction: shearlets, radial 



(f) Reconstruction: wavelets, radial 



(li) Reconstruction: Fourier inversion, radial 


Figure 4: Reconstruction using a spiral and radial mask. 
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7 Future work 


We like to address multiple extensions and future directions of the work presented in this paper. 

7.1 Linear stable sampling rate 

In this paper we have established a stable sampling rate for sampling with respect to complex 
exponentials and recovery using compactly supported shear let frames. Our rate essentially dif¬ 
fers from the rate for wavelets by an additional N 5 term and the lower frame bound. Although 
it seems plausible that the rate should depend on the stability of the redundant system, mea¬ 
sured by the lower frame bound, it is interesting to see, whetherthe could be improved to actual 
linearity up to a dependency of the lower frame bound, i.e. erasing N s while keeping An- 
Furthermore, efficient numerical algorithms to compute the stable sampling rate for shearlets 
are currently missing. This issue is left for future work. 

7.2 Nonuniform sampling 

Moreover, we have restricted ourselves to the case of a uniform sampling pattern, in fact, the 
samples were assumed to be acquired on a grid. Although this is possible in practice, it is 
very unusual. Many MRI acquisition schemes are nowadays non-uniform schemes, e.g. spiral 
patterns or radial patterns are often used, see Figure 0a] and Figure I4bl In [2] [Tj the authors 
introduced non-uniform generalized sampling and proved stable reconstructions are possible 
when using folded, periodic, and boundary wavelets. We strongly believe that the methods used 
in the present paper can be used to obtain similar results for shearlets in the non-uniform case 
as well. This is will be studied in an upcoming work. 

7.3 Parabolic molecules 

Shearlets are a special type of parabolic molecules 0G2- It is therefore an interesting question 
whether our results can be extended to such systems. We explain briefly how such generaliza¬ 
tions can be tackled. 

First of all, the reconstruction space 77/v is spanned by the first N elements of the repre¬ 
sentation system. In our case of compactly supported shearlets it was build by considering all 
atoms whose support intersect a region of interest. If one considers atoms that are not com¬ 
pactly supported, then one has to make a careful choice of ordering the representation system 
and select the first N elements of the system. However, the argument for when a reconstruction 
from Fourier measurements is possible, is solely based on very good frequency localization of 
the atoms. Thus this argument is applicable to other systems such as wavelet frames or general 
parabolic molecules. However, for a qualitative judgement of the results one needs to relate the 
lower frame bound of the finite system to the approximation rate, similar as we did in (14.4|) . 

7.4 Compressed sensing 

Another perspective that is indispensable for reconstruction problems from Fourier measure¬ 
ments such as in MRI is compressed sensing cf. [HB El- Moreover, the new concepts of 
asymptotic sparsity and asymptotic incoherence introduced in [8] are perfectly suited for shear- 
lets, in fact natural images are asymptotically sparse in shearlets and it can be shown that 
the cross gramian between complex exponentials and shearlets is asymptotically incoherent. 
Further investigations in this direction are also possivle for future work. 

As shown in the numerics, shearlets have a great potential in recovering real medical im¬ 
ages from their subsampled Fourier data. The development of a sophisticated l 1 algorithm is 
therefore of high importance. A mathematical analysis and a comparison of different recovery 
algorithms will be provided in an upcoming paper. 
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